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SUMMARY 

Solutions to the governing equations of a spur gear transmission model, 
developed in NASA TM-100180 (AVSCOM TM-87-C-2) are presented. Factors affect- 
ing the dynamic load are identified. It is found that the dynamic load 
increases with operating speed up to a system natural frequency. At operating 
speeds beyond the natural frequency the dynamic load decreases dramatically. 
Also, it is found that the transmitted load and shaft inertia have little 
effect upon the total dynamic load. Damping and friction decrease the dynamic 
load. Finally, tooth stiffness has a significant effect upon dynamic loading: 
the higher the stiffness, the lower the dynamic loading. Also, the higher the 
stiffness the higher the rotating speed required for peak dynamic response. 


INTRODUCTION 

The development of a simple parallel shaft spur gear transmission model 
with its dynamic differential equations and solution procedures were presented 
in reference 1. Various parameters such as inertia, stiffness, friction and 
damping were included in the governing equations for further study (ref. 1). 

The purpose of this report is to determine the effect of these parameters 
on gear dynamic load. The dynamic load of a gear transmission can be found by 
solving the governing equation. The solution is known as the dynamic motion 
of the gear transmission. This dynamic motion can then be substituted into 
other analytical formula and solved for gear dynamic loads. 


‘Member, ASME. 


A model of the transmission is depicted in figure i. 
are: 

,3 M®M + C sl <0 M * V + K s 1 <0 M " V = T M 


The governing equations 

(1) 


Ji e i + C s1 <e r e M > + k $ 1 (0 1 - e M ) + c g (t) [R b1 e - R b2 e 2 ] 
+ K g (t) CR b1 (R bl e r R b2 e 2 )] = T f1 (t) 


( 2 ) 


J 2 0 2 + C s2 (0 2 0 1 } + K s 2 <0 2 “ 0 1 > + C g (t) [R b2 0 2 " R bl 0 l ] 

+ K g (t) [R b2 (R b2 e 2 - R b1 e ] ) ] = T f2 <t) 


(3) 


J L 0 L + C s2 <0 L " 0 2 ) + K s2 (0 L " V = “ T L 


(4) 


where J-j , J 2 , and J|_ represent the mass moments of inertia of the motor, 
the gears, and the load; C s ] , C s2 , and Cg ( t ) are damping coefficients of the 
shafts and the gears; K s ], K s2 , and K g (t) are stiffness of the shafts and the 
gears; T^, T|_ , Tf](t), and Tf 2 (t) are motor and load torques and frictional 
torques on the gears; R b i and R b2 are base circle radii of the gears; t is 
time, and the overdots indicate time differentiation. 

In this report we present the results of numerical solutions of equa- 
tions (1) to (4) for a typical transmission system. A flow chart outlining 
the numerical procedure is presented. Natural frequencies are determined. 

The dynamic load is determined. Finally, the dynamic factor, defined as the 
ratio of the dynamic load to the static load, is determined.' The results are 
calculated as functions of the rotating speed and roll angle for a variety of 
damping and stiffness conditions. 


NOMENCLATURE 

Cg damping coefficient of gear tooth mesh, N-sec (lb-sec) 

C s damping coefficient of shaft, N-m-sec (in. -lb-sec) 

J|_ polar moment of inertia of load, kg-m 2 (in. -lb-sec 2 ) 

Jm polar moment of inertia of motor, kg-m 2 (in. -lb-sec 2 ) 

0] polar moment of inertia of gear 1, kg-m 2 (in. -lb-sec 2 ) 


Uhe term "dynamic factor" or "dynamic load factor" has been used inconsis- 
tently in the literature. The American Gear Manufacturer's Association 
(AGMA) dynamic factor K v is used as a strength reduction factor and is 
defined as the maximum static load divided by the maximum dynamic load. This 
paper will follow the ISO convention which uses the dynamic factor as a 
load/stress increasing factor. Therefore, Kj = 1/K V . 
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J2 polar moment of inertia of gear 2, kg-m 2 (in.ib-sec 2 ) 
K(j dynamic factor 

Kg stiffness of gear tooth, N-m/rad (in. -lb/rad) 

K s stiffness of shaft, N-m/rad (in. -lb/rad) 

K v AGMA dynamic factor, Kv = 1 /K<j 
R5 base radius, mm (in.) 

Rp pitch radius, mm (in.) 

T|_ torque on load, N-m (in. -lb) 

torque on motor, N-m (in. -lb) 

Tfl torque on gear 1 , N-m (in. -lb) 

Tf2 torque on gear 2, N-m (in. -lb) 

W applied load, N (lb) 

0 angular displacement, rad 

0 angular velocity, rad/sec 

0 angular acceleration, rad/sec 2 

u n natural frequency, Hz 

E damping ratio 


PROCEDURE 

Figure 2 presents a flow chart of the computational procedure used in the 
parameter study. The procedure is the same as that outlined at the end of 
reference 1 . 

In conducting the analysis it is useful to compare the locations of the 
peak dynamic loads with the locations of the system natural frequencies (or 
critical speeds). The natural frequencies themselves may be obtained by exam- 
ining the undamped equations of motion. These equations may be written in the 
matrix form: 


[J] [ 0 ] + [K] [ 0 ] = [ 0 ] 
where the inertia matrix [J] is 

0 0 0 

J 1 0 0 

0 J 2 0 

0 0 j l 



( 5 ) 


( 6 ) 
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and the stiffness matrix C K] is 


CK] = 


si 

- K sl 


si 

<yavg R2 bl 

-<v 


-(K g ) avg R bl R b2 

K sl * 


Vavg R b 2 

- K s2 


0 

0 

-K 


s2 


's2 


( 7 ) 


where (kg) aV g represents the average value of the gear mesh stiffness. It 
is taken as tne sum of the discrete tooth stiffness values of a mesh cycle 
divided by the number of mesh positions in the cycle. 


In the parameter study the system had identical gears with the following 
properties : 


Number of teeth 

Module, mm 

Pitch diameter, m (in.) 

Pressure angle, deg 

Applied load, N (lb) 

Face width, m (in.) 

Moment of inertia, kg m 2 (in. -lb sec 2 ) 

Average tooth stiffness, N m/rad (lb-in. /rad) . . . 
Damping ratio 


36 

3.18 (8 diametral pitch) 

0.1143 (4.5) 

20 

2670 (600) 

0.0254 (1.0) 

. 3.3323x10-3 (0.02947) 

3.991x10 s (3.5355x10 s ) 
0.10 


The shaft stiffness and inertias were: 


Shaft stiffness, N m/rad (in. -lb/rad) 1138.17 (10081) 

Load inertia and motor inertia, 

kg m 2 (in. -lb sec2) 9.989x10-3 (each) (0.08841) 


RESULTS 


Using the aforementioned data in the gear system model shown in figure 1, 
the natural frequencies of the four degrees of freedom model were found to be 


o n1 = 0 (rigid body mode) u n2 = 1.49 Hz 

0 ) n3 = 2.99 Hz co n4 = 144.8 Hz 


( 8 ) 


The first three natural frequencies are well below tooth meshing frequencies 
and are therefore not of interest in this analysis (although they can still be 
excited by shaft eccentricity which is not modeled here). The fourth natural 
frequency matches tooth meshing frequency at the critical speed of 8688 rpm 
which is within the operating range of the gears. 


Figure 3 shows the variation of dynamic load response for a pair of teeth 
as a function of roll angle. At speeds much lower than the critical speed, 
the dynamic load response is basically a static load sharing in phase with the 
stiffness change, superimposed with an oscillatory load at a frequency corre- 
sponding to the natural frequency. 
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At higher speed, close to the critical speed, the dynamic load variation 
becomes so abrupt that it produces tooth separation. The peak dynamic load is 
much higher than the static load and is very likely to be a source of gear 
noise and early surface fatigue. 

Figure 4 shows the dynamic factor K<j as a function of operating speed. 
Prominent peaks (resonances) may be seen at speeds of 7650 and 4200 rpm. The 
larger peak 7650 rpm occurs at 88 percent of the calculated critical speed. 

The experimental work by Kubo (ref. 6) reported a similar result that the crit- 
ical speed was found at about 90 percent of the calculated critical speed. 

The second peak at 4200 rpm is a nonlinear effect of the time varying tooth 
stiffness known as the parametric resonance. This parametric resonance fre- 
quency occurs at about one-half of the critical speed (ref. 2). For speeds 
above the critical speed, the dynamic response decreases steadily in the same 
manner as with elementary vibrating systems. 

Figure 5 shows a three-dimensional representation of the system dynamic 
response. The horizontal axis represents the operating speed, and the contact 
position along the tooth profile. The total number of contact positions is 
121. The vertical axis is the dynamic factor. 

Figure 6 presents a contour plot of the system dynamic response. The 
shaded areas represent regions where tooth separation occurs. They are locat- 
ed in the double contact regions. At near resonance speeds the vibration 
amplitude exceeds the deflection of the meshing teeth, thus inducing tooth 
separation. 

As the speed increases, the dynamic response also shows a phase shift 
towards the higher numbered contact positions. This phenomenon can be seen by 
noting that at speeds from 600 to 4200 rpm (one-half subharmonic), the maximum 
dynamic load occurs at double- to single-contact transitions (position 51) and 
gradually changes to single- to double-contact transitions (position 75). At 
speeds between one-half subharmonic and resonance the maximum peak stays near 
the single- to double-contact transitions. After the speed passes resonance, 
the major dynamic peak moves again towards higher contact positions on tooth 
profile with increasing speed. 

Figure 7 shows the dynamic factor as a function of the transmitted loads 
for three different speeds. There is only a small decrease in dynamic factor 
with increased applied load. 

Figure 8 shows the effect of damping on the dynamic load. It is seen 
that damping has its greatest effect near resonance frequencies. 

Changes in shaft stiffness have a minor effect on the system dynamic 
response. However changes of tooth stiffness have a major effect on the 
response. Figure 9 shows that the higher the tooth stiffness the lower the 
dynamic response (dynamic factor). This is consistent with observations that 
as the tooth stiffness increases the effect of gear mass on the system dynam- 
ics is reduced. Figure 9 also shows that system resonance frequencies are 
increased as the tooth stiffness increases. This is a potentially useful 
effect for the design of gear systems. 

The effect of shaft mass is assumed small compared to that of the gears. 
Figure 10 shows that as the gear inertia is reduced, the dynamic response is 
also reduced. 
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For gears with different diametral pitches, the dynamic response is dif- 
ferent due to the change in contact ratio. Gears with a finer pitch have a 
higher contact ratio. Since the contact ratio is a measure of the duration of 
the load being shared by more than one pair of teeth, it has a significant 
effect on the system dynamic response. 

Figure 11 shows a comparison between gears having different diametral 
pitches. The finer pitch gears, having a higher contact ratio, have a smaller 
dynamic load than the coarser pitch gears. 


DISCUSSION 


In 1927, A. A. Ross (ref. 
for the dynamic factor k v : 


3) introduced the following empirical formula 


k v = 


78 


(78 + 


yv) 


(9) 


where v is the pitch line speed measured in ft/min. This expression received 
acceptance as a standard factor used by the American Gear Manufacturer's 
Association (AGMA). In 1959, a similar factor for use with higher precision 
gears was introduced by Wellauer (ref. 4) 


k = 78 

V 78 + yv 


( 10 ) 


Equations (9) and (10) are recognized as being conservative when applied 
with very high precision gears. They are thought to predict dynamic loads 
which are higher than the actual loads. 

Buckingham (ref. 5) has also developed an expression for the dynamic load 
in terms of the pitch line speed and the applied load. His formula is 

«d ■ W * [f a< 2 V V 1 * <"> 

where is the dynamic load, W is the applied load and the factors f a 
and fb are 

f = f.f /(f, + f ) and f. = 0.0000555 EF + W (12) 

a b c b c b 

where F is the face width, E is the elastic constant, and f c is 

f = 0.00120 C(R 1 + R 2 )/R ] R 2 ]mv 2 (13) 

where Ri and R 2 are the pitch radii of the gears and m is their effec- 
tive mass. (In these expressions the units are in pounds and inches except 
for the pitch line speed which is measured in ft/min.) 
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Kubo (ref. 6) measured static and dynamic gear stresses on several 
high-precision spur gear systems. Kubo expressed the dynamic factor as the 
ratio of maximum dynamic to the maximum static stress. Since stress is propor- 
tional to load, Kubo's definition of dynamic factor is identical to that used 
here. Figure 12 shows a comparison of (1/K V ) from the AGMA high-precision for- 
mula (eq. (10)), (W(j/W) from Buckingham's formula, Kubo's results, and the 
results of the computer simulation for an identical spur gear pair with 4 mm 
module, 25 teeth, 20° pressure angle, 15 mm face width, 131.5 kN/m applied 
load, and 207 GPa Young's modulus. 

There is a good agreement between Kubo's result and the computer 
simulation. 


CONCLUSIONS 

From the foregoing results several conclusions may be stated: 

1. For accurately made gears, the dynamic load is significantly affected 

by the contact ratio: for increased diametral pitch, that is, for high contact 

ratio gears, the dynamic load is decreased. 

2. The tooth stiffness has a significant effect upon the dynamic load: 
the higher the stiffness the lower the dynamic load. Also, the higher the 
stiffness the higher the critical speed for peak response. 

3. The dynamic load generally increases with the operating speed until a 
resonance is reached. The dynamic load decreases rapidly beyond the resonance. 

4. Damping and friction decrease the dynamic load with the most dramatic 
effects occurring near the resonance frequencies. 

5. The applied load has a relatively minor effect upon the dynamic factor. 

6. Tooth separation — leading to impact — occurs in the double tooth 
contact region since the deflections are smallest in that region. 

7. The dynamic factor is largest for contact points near the tooth tip. 

8. The dynamic factor decreases with decreases in the gear body inertia. 
Shaft moment of inertia has a minimal effect upon the dynamic factor. 
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FIGURE 1. - A SIMPLE SPUR GEAR SYSTEM. 


GEOMETRY OF SYSTEM COMPONENTS AND 
CONDITIONS OF SYSTEM OPERATION 



OUTPUT OF RESULTS 


FIGURE 2. - FLOW CHART OF COMPUTATIONAL PROCEDURE. 
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FIGURE 3. - DYNAMIC LOADS AT DIFFERENT SPEEDS FOR IDENTICAL YEARS. FIGURE 4. - DYNAMIC FACTOR AS A FUNCTION OF ROTATING SPEED. 

MODULE = 3.18 MM; PRESSURE ANGLE = 20°; PITCH RADII = 57.1 MM; 

APPLIED LOAD = 2760 N (600 LB). 


2.5 r 



2200 

ROLL ANGLE INDEX SPEED, RPM 
FIGURE 5. - SYSTEM DYNAMIC RESPONSE. 
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FIGURE 7. - EFFECT OF APPLIED LOAD ON DYNAMIC 
RESPONSE. 
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FIGURE 8. - EFFECT OF DAMPING ON DYNAMIC RESPONSE 
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FIGURE 10. - EFFECT OF GEAR INERTIA ON DYNAMIC ON DYNAMIC FIGURE 11. - EFFECT OF DIAMETRAL PITCH ON DYNAMIC RESPONSE. 

RESPONSE. (GEAR INERTIA: 100% VALUE = 3.33 X 10" 3 KG-M 2 ). 
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NORMALIZED SPEED, w/U n 

FIGURE 12. - COMPARISON OF EXPERIMENTAL, EMPIRICAL, AND 
COMPUTER SIMULATION RESULTS. 
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